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Abstract 

We develop a new method to calculate eigenvalues in frustrated quantum 
spin models. It is based on the stochastic state selection (SSS) method, which 
is an unconventional Monte Carlo technique we have investigated in recent years. 
We observe that a kind of equilibrium is realized under some conditions when 
we repeatedly operate a Hamiltonian and a random choice operator, which is 
defined by stochastic variables in the SSS method, to a trial state. In this 
equilibrium, which we call the SSS equilibrium, we can evaluate the lowest 
eigenvalue of the Hamiltonian using the statistical average of the normalization 
factor of the generated state. 

The SSS equilibrium itself has been already observed in unfrustrated mod- 
els. Our study in this paper shows that we can also see the equilibrium in 
frustrated models, with some restriction on values of a parameter introduced 
in the SSS method. As a concrete example, we employ the spin-1/2 frustrated 
J1-J2 Heisenberg model on the square lattice. We present numerical results on 
the 20-, 32-, 36-site systems, which demonstrate that statistical averages of the 
normalization factors reproduce the known exact eigenvalue in good precision. 
Finally we apply the method to the 40-site system. Then we obtain the value 
of the lowest energy eigenvalue with an error less than 0.2%. 

1 Introduction 

The quantum Monte Carlo method is well established in the numerical study of quan- 
tum spin systems. This method has given us fruitful results to understand properties of 
unfrustrated systems, especially of the spin-1/2 quantum Heisenberg anti-ferromagnet on 
bipartite lattices [1]. But, due to the so-called sign problem, the method is not useful 
for studies of frustrated systems. In contrast to this the exact diagonalization is applica- 
ble even when systems are frustrated. By this method, however, one can only deal with 
systems on small-sized lattices. In order to overcome this difficulty lots of active studies 
for numerical methods have been made in the Monte Carlo approach. Among them the 
re-normalization group method [2] and the reconfiguration method combined with the 
fixed node method P] are quite noticeable. Also, the extensions of the density matrix 
re-normalization group method |1] is worth noting. 
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Recently we have developed a new Monte Carlo method to evaluate energy eigenvalues 
of quantum spin systems El El E]- We call it the stochastic state selection (SSS) 
method |llO|. This method provides an essential improvement on early attempts to trun- 
cate the Hilbert space by omitting small values of coefficients [11] . Not being based on 
importance samplings, the method enables us to select a relatively small number of states 
from a vast full vector space. This selection is simple and is mathematically justified so 
that one can calculate correct values of inner products through statistical averaging pro- 
cesses. In order to realize this selection, we employ an operator which we call the random 
choice operator. This operator is represented by a diagonal matrix whose elements are 
stochastic variables with unit averages. 

It is possible to combine the SSS method with various techniques for numerical studies. 
For instance, several applications with the power method are in refs. |6l [7] . By repeating 
alternate operations of the random choice operator and the Hamiltonian to a trial state, 
we numerically demonstrated that we can obtain expectation values for powers of the 
Hamiltonian even when limited computer memory resources are available. We also have 
combined the SSS method with the Lanczos method in order to obtain the lowest energy 
eigenvalue of the spin-1/2 quantum Heisenberg anti-ferromagnet on the triangular lattice 
up to 48 sites [9J. It should be noted that these applications of the SSS method are simple 
from mathematical points of view, because the statistical averaging process for random 
variables is clear. 

In this paper we discuss another application of the SSS method that is essentially dif- 
ferent from the applications stated above. Here we consider a number of intermediate 
states which are successively generated by the random choice operator and an operator 
related to the Hamiltonian, and measure their normalization factors. Generally the stan- 
dard deviation of the inner product between the initial trial state and the L-th generated 
state increases as L grows. In some conditions, however, the deviation does not increase 
but become constant when L goes beyond some value. We call this phenomenon as an 
equilibrium, because in this situation statistical averages of the inner product do not de- 
pend on L after suitable normalizations. This type of equilibrium is quite interesting 
since from values of the normalization factor one can obtain the lowest energy eigenvalue 
precisely and simply. We first observed such phenomena in our study of the J1-J2 Heisen- 
berg model with unfrustrated [12] couplings [8]. In this paper we study the equilibrium 
in the frustrated case. From qualitative discussions we show that the equilibrium for 
frustrated systems exists in some parameter region. A concrete example for numerical 
examinations is the spin-1/2 J1-J2 Heisenberg model on the square lattice with couplings 
Ji = 1 and J2 = 0.5, the model which is one of the most popular frustrated models with 
ample numerical results [131 El ESI EH Ell- Studying this model up to 40 sites, we find 
that our results give tangible evidences for the equilibrium. 

This paper is organized as follows. Next two sections are devoted to definitions and 
qualitative discussions. In Section [2l after brief descriptions for the SSS method, a defi- 
nition of the SSS equilibrium is given. Here we also explain why the equilibrium enables 
us to calculate the lowest energy eigenvalue. Section [3] is for summarized discussions on 
where the equilibrium takes place. We briefiy repeat our previous argument [8J with which 
we concluded that the equilibrium should appear in any unfrustrated system. Then we 
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argue that for frustrated systems the equihbrium also exists in some parameter region. 
In Section H] we present numerical results on the frustrated J1-J2 Heisenberg model. We 
first direct our attention to the model on a 20-site lattice. On this small lattice it is 
easy to obtain the exact ground state by the diagonalization method. Consequently we 
can directly measure overlaps between the exact ground state and the states generated 
through stochastic state selections. Note that, as is discussed in Section |2l those overlaps 
are most effective to prove the existence of the SSS equilibrium. In the 20-site case it 
is also possible to make detailed observations on several quantities introduced in Section 
[31 In Subsection 14.11 we will show that assumptions we employ in Section [3] are mostly 
reasonable for the model and that the eigenvalue we obtain through normalization factors 
is in good agreement with the exact ground-state energy. We also examine how many 
basis states are needed in our method to calculate the eigenvalue with enough precision. 
Subsection 14.21 is to present our results for the 32- and 36-site systems. For these sizes we 
construct bases taking some symmetries into account. Without knowing the exact ground 
state, we use not overlaps but normalization factors to observe the SSS equilibrium in 
these systems. Our results indicate that the eigenvalues we calculate in the equilibrium 
nicely reproduce the exact ground-state energies reported in ref . [13] . In Subsection 14.31 
we present results for the 40-site system. The final section is devoted to summary and 
discussions. 



2 Stochastic State Selection Equilibrium 

First we briefly review the SSS method [7j. Consider a system whose Hamiltonian is H 
and suppose its full vector space is given by a basis {| i)} {i = 1, A^)- Here iVy denotes 
the size of the full vector space. We denote the lowest energy eigenvalue of if by -E and 
its eigenstate by | iPe)- For convenience a new notation Q is used hereafter, 

Q = li-H, (1) 

with the identity operator / and a positive number I. Since the present method is based 
on the power method, the value of / should be chosen to ensure that Q = I — E{> 0) is 
the largest among absolute values of eigenvalues for Q. 

Let us employ a normalized trial state | tp^^'') = J2i I i)cf'\ (V'''^'* I 'ip^'^^) = 1- We 
introduce random choice operators 

M{,('")} = EM)#^(^I (m = l,2,---), (2) 

i=l 

in order to calculate 

QM{^(L)}QM|^(L-i)} ■ ■ ■ QM{^(i)} I V-^"^) (3) 

instead of \ ip^'^^), the state which would become close to (/ — E)^ \ ipEj for large 
values of L. The random variable r/J""^ in ^ is generated following the on-off prob- 
ability function P^^\r]). Using a positive parameter e and the coefficient c|™~^'* in 
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the normalized intermediate state | t/'*-™' ^'') = I ^c-™ ^\ which is proportional to 

QM^^(m-i)j ■ ■ ■QM^^(i)y I ip^'^^), we define 



P. 
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(m) 
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(4) 
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Then starting from a given | il)'^'^'>) we sequentially calculate 

I ^(™)) = QM^^M} I (5) 
for m = 1, 2, 3, ■ ■ ■, where C^"^-* (> 0) is the normalization factor determined by 

= I M{^(™)}Q2m^^(™)} I (6) 

We also define a random state \ with the normalization condition (x^™'* | X^™'*) = 

Iby 

i 

It should be kept in mind that for any state | = | i)hi with 6j's irrelevant to {fif^^^^} 
the statistical average of the inner product between | $) and this random state is zero, 



because 



from (jl]) [IBJ. Also remember that 
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where we use (El) and 
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Now we define the SSS equilibrium [19j. We divide the intermediate state | ip^^'>) into 
a part which is proportional to | V'e) and the rest, 



Here 



(12) 
(13) 
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and 



with the normahzation condition (C*-™'' | C*'™'') = 1- Note that 
and 



w 



(m) 



+ 



(m) 



(14) 
(15) 
(16) 



by definition. What we mean by the SSS equihbrium is that there exists a hmit w'^^'^^ 
defined by 



ms+mt-l 

hm — 

TTlt^OO 777, ^ ' 



(0 < < 1 



(17) 



where w^'^'^^ is independent of rUs whenever nis is greater than or equal to some value of 
m. 

Finally let us comment on a relation useful in the equilibrium in order to extract the 
value of E = I — Q. Using ([5]), (I7j) and ( fT3i) we obtain 

_ I ^(m+1) 

= (^E I Q { I I } 

= Q{ (^E I ^('")) + (^E I X^™^)^?^'"^ } 

= g^('") + Q(^E|x(™))^?(™^ (18) 

If the second term in the right-hand side is negligible, it leads to 



w 



(m+l) 



Q 



(m) 



C(m+1) ■ (19) 

Then, in the equilibrium where w^™) ^ ~ y;(eq) ^ g fQj^ sufficiently large values of 

m, we can expect 

Q ^ c'^rn+1) _ ^20) 

Thus we become aware that the value of E can be estimated from the normalization factor 



3 Existence of the SSS Equilibrium 

In this section we present analytic and qualitative discussions on the existence of the 
SSS equilibrium. This section consists of three subsections. The first subsection gives 
some equations for C^"^~^^^ assuming that fiuctuations can be neglected. Next subsection 
is devoted to summarized discussions for unfrustrated systems, where all elements of the 
operator Q are non- negative [5]. In the third subsection we show that the SSS equilibrium 
also exists for frustrated systems when the parameter e is small enough. 
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3.1 An Equation for the Equilibrium 

We pay our attention to a relation led from ([5]), ([7j) and (fT2!) . 



(21) 



With normalization conditions and (fT5l) and (fT6|) it yields 



^(m+l) 



w 



(m) 



(m)\ 



(m) 



(m)\ 



(m) 



(22) 



We assume that both | | (^(™')) and (x''™'' I I X*^™^) ^ire independent of m when 
m is large enough, 



Q2C 



(23) 
(24) 



where ^-iid denote positive constants. It should be noted that both Q2C Q2X 
are always less than Q"^ because = {ipE I I V^e) is the largest one among {ijj \ Q"^ \ ip) 
by definition. We also assume both of cross terms in fl22p are negligible, notifying that 
statistical averages of them vanish. Thus we obtain a relation 



^(m+l) 



Let us here emphasize that 



w 



(m) 



(m) 



+ Q2C 1 - 



W 



(m) 



Q2x 



-,{m) 



(25) 



measures the degree of the difference between | tp^"^^) 



and M. 



{^(m + l) } 



ipi"^)'^. We would like to notify again that the statistical average of 



(m) 



is calculated from coefficients for | -j/;*^™)) and the parameter e (see f llUj) ). 
3.2 The Equilibrium in Unfrustrated Systems 

Here we brieffy comment how we conclude that the SSS equilibrium should exist in unfrus- 
trated cases with any value of e. Detailed discussions and several numerical examinations 
with the J1-J2 Heisenberg model {J2/J1 = 0, —1) are presented in ref |8]. In this sub- 
section we limit ourselves to the case | ip^^^) =\ ip^) in order to make our analysis clear. 
Discussions for a good approximate initial trial state can be made in a similar manner. 

For unfrustrated systems it is always possible to choose an adequate basis {\ i)} for 
which all /j's in | ip^) = J2i I i)fi as well as qij = {i \ Q \ j)'s are non-negative. Note that 
all coefficients in the expansion of | ip^^^), \ '4'^'^^)^ " " " are also non-negative then. This is 
because, in the relation we learn from 



c 



(m) 



(26) 
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Qij > for all i and j and rjj^^ /C^"^'' > for all j by definition. Let us then examine the 
first term in (ITU]) . One upper bound for it is given by 



e E l^'"^! < ^El4™^l=^E4"^^ (27) 

where the last equality follows from the fact that all c\"^^ are non-negative here. Further, 
we can expect that, in positive definite cases, 

eEc^^-^-^'-^eE/^ (28) 

i i 

holds for any m [8j. Now we add one more assumption in order to clearly propose a 



relation to be expected in the SSS equilibrium. Noting that Y^r 
for the second term of ffTOl) that 



1, we assume 



E [4"^T ^ ^ (29) 



\c\ '\<e 



holds with sufficiently large values of m, where i^' (0 < i^T < 1) is a constant defined by 
e. This brings, together with fl28|) . a relation 



Gw^^'^-K, G = eY^U. (30) 



Combining flTI?]) . fl2^ and fl5Ul) . we obtain a recursive relation for ty^™"), 

. (31) 

V(Q' - Q2c) + GQax^^"^) + (Q2C - ^Q2x) 

A simple mathematical analysis on fl3T|) |8] leads us to the conclusion that w^'^'^^ exists 
as far as G > i^' and (52C > ^Q2x 120] ■ Also it is concluded that w^'^'^^ should satisfy a 
quadratic equation for w, 

(Q^ - Q2c) + GQax^ - (Q' - Q2C + KQ^x) = . (32) 
The relevant solution which belongs to the interval (0, 1] is 



Thus we come to a conclusion that the SSS equilibrium exists in unfrustrated systems. 
3.3 The Equilibrium in Frustrated Systems 

In this subsection we show that the SSS equilibrium should, at least for small values of e, 
also exist for frustrated systems. 
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Let us find an equation for w^'^'^'' in frustrated systems from (l25l) . We again start 



from examining 



(m) 



with an assumption that it becomes almost constant (= g"^) for 



sufficiently large values of m. Then (fTOl) brings an upper bound, 



g^ ~ 



(m) 



(m) 



')) < ^ E 

I 1 



.Ml 



< ^E 



.Ml 



(34) 



< e 



i=l 



The last inequality is based on the fact that max E l*^*! " 

If the SSS equilibrium exists, w^'^'^^ should satisfy the following equation for 



{Q' - Q2c) (i 



w 



(35) 



which we obtain using (I2U]) and (1251) . Clearly this equation, where — Q2c_ > and 
Q2x — always hold by definition, has a relevant solution if Q2x9^ is less than — Q2C- 
Because of the upper bound flMl) we then see that there exists a solution 



w 



(eq) 



_ Q2x9^ 



(0 < < 1) 



(36) 



for small values of e. Therefore we conclude that the SSS equilibrium is realized in 
frustrated systems when the value of the parameter e is small enough. 



4 Numerical Study 

Now we present numerical results on the frustrated J1-J2 Heisenberg model. The Hamil- 
tonian of the model is 

HjiJ2 = Ji'^ Si ■ Sj + J2 E ' ^3' ■ ^'^^) 

{nn) {nnn) 

Here Sk denotes the spin 1/2 operator on the site k and summations run over the nearest 
neighbor pairs (nn) or over the next-nearest neighbor pairs {nnn) of the square lattice. 
In this work we restrict ourselves to the = sector, where 5*2 denotes the z component 
of the total spin [21]. Values of couplings are fixed to be Ji = 1 and J2 = 0.5 throughout 
this paper. 

As is mentioned in Section O the value of / in the operator Q = 11 — Hj^j^ should be 
determined to make the largest eigenvalue for Q correspond to the lowest eigenvalue of 
Hjij2, which we denote by E. Since the largest eigenvalue of Hj^^j^ on an A^^-site lattice is 
2( Ji + J2)Ns/4: with the largest Sz{= Ns/2), this means a condition 21 > ( Ji + J2)Ns/2 - 
\E\. Also keep in mind that very large value of / should be avoided because it causes slow 
convergence. Values of / thus chosen for each lattice size with a guess for E will be given 
in each subsection. 

This section, which includes three subsections, is for numerical studies on lattices up 
to Ns = 40. In the first subsection we make a detailed study of the SSS equilibrium on 
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a. Ng = 20 lattice. Comparing our results with those obtained from the exact eigenstate, 
we demonstrate that our assumptions described in the previous section are reasonable. 
In the second subsection we evaluate the lowest energy eigenvalue on Ng = 32 and 36 
lattices, imposing some symmetries on their bases. We will see that our results are in 
good agreement with the known exact ground-state energy. The third subsection is for 
the Ns = 40 results. 



4.1 Ns=20 Results 

For this lattice size there are 2oCio=184,756 configurations which fulfill the condition 
Sz = 0. We use these configurations as basis states | i), without assuming any symmetry 
on the state | iIj^"^^). Let Q = 31 — Hj^j^ in this subsection. It is easy to carry out 
the exact diagonalization of the matrix [(i | Q | j)] and find the exact eigenstate | t/'e)- 
Remember that Q \ i/je) = Q I ^e) and Q is the eigenvalue for which the value of \Q\ 
is the largest among all eigenvalues for Q. Numerically we obtain Q = 13.0123, which 
certainly reproduces the exact ground-state energy of the system. 

In this subsection we start from the exact eigenstate, namely, | ip^^^) =\ -ip^). Then 
according to ([H]) we sequentially calculate | -ip^"^^) up to m = 200. Figs. [T] - [H] present 
results on several quantities described in previous sections. All statistical averages in 
these Figures are calculated from 10 samples for each value of e. We do not plot statistical 
errors because they are scarcely beyond the marks we used in the Figures. 

In Fig. [1] we present the results on w^"^^ = {x/je \ ip^"^^) for several values of e as a function 
of m. We see that, for each value of e less than or equal to 0.02, w^"^^ {m> 30) fluctuates 
around a finite value defined by e. This is the very evidence for the SSS equilibrium. 
While, when e > 0.025, w^*") rapidly decreases to zero as m grows. Therefore no SSS 
equilibrium takes place in this parameter region. Fig. [2] shows the normalization factor 
(ji'fn) calculate from (jS]). Our results in Fig. [2] confirm that values of C*^™-' are around 
the exact eigenvalue Q when the SSS equilibrium occurs. 

Figs. [31 m and [5] are to endorse relations fl23|) . flM]) and fl25l) . In Fig. [3] we plot values 
of ((^(™') I Q"^ I C*^™)) and (x*-'"'* I I X*'™''') foi' some values of e. It is clear that our 
assumptions ([23]), ([21D are both acceptable because (C^™^ | I C^"^) and (x^""^ | | 
X^^^) hardly depend on m when m > 20. We also observe that fl23|) and fl24|) hold 
for other values of e. The resultant value of {Q2x) fi'^st increases (decreases), then 
decreases (increases), as we lessen the value of e. The change for is, however, too 
small to read in this Figure. Figs. [H and [5] present ratios of three terms in the right-hand 

side of fl25p to (7^™+^-' for e = 0.025, where no equilibrium is observed, and for e = 0.02 
with which the system realizes the equilibrium. We see that the sum of the three ratios 
for each m is in good agreement with the expected value, namely 1.0. We also observe 
such agreements for other values of e. Therefore we can say that the relation fl25|) is 
acceptable regardless of the equilibrium. In other words, we can conclude that the cross 
terms 2Q^w'■'^^^pE \ x(™))c/(™) and 2s'■"'^C^"''^ \ \ x^""'^) g'-'^^ in ([22]) are always negligible 
for this system. 

r -|2 

Let us next examine the relation f[36]) in the SSS equilibrium. Fig. [6] plots (( g^'"^' )) 
with e = 0.02, 0.015, 0.01 and 0.005. We see, in accordance with what we have assumed 
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at the beginning of that fluctuations in g'^"^^ due to different random number se- 

r . s-|2 

quences are neghgibly small and that g^"''' for each value of e becomes almost irrelevant 
to m for m > 20. For any value of e shown in Fig. [6] we observe that the value of w^'^'^^ 

estimated from ^1 - (( (x^"") I I X^"^) )) (([c/^™^]^))/ (Q^ - (( (C^""^ I I C^""^) ))) is in 

good agreement with values of ((w*^™-')) in Fig. [TJ 

Finally we comment on how many basis states should be included in our calculation. As 
a result of the stochastic state selection, A^,^™-*, the number of non-zero coefficients in the 
expansion of M|^(m+i)} | '?/'*^™^), is much less than A^b™"*, the number of non-zero coefficients 

in the expansion of | tjj'^"^^). Hence it is not A^^™-* but A^™"* that determines the necessary 
memory resources. For each value of e we observe that after several iterative operations 
to generate QM{^(™)} | i/j''"^~^'^){= C'^"^^ \ ?/'^™))), A",^"*^ does not increase anymore. The 
data also show very small deviations for different random number sequences. Therefore 
the upper limit of nI^"^\ which is dependent on the value of e, is critical for the numerical 
study. Some results for the 20-site system are as follows. When e = 0.05 the upper 
limit of A^,^'"'* is slightly less than 1 x 10^, namely we can do with only one half of the 
A^v = 184, 756 basis states. With the value e = 0.02, for which the system exhibits the 
SSS equilibrium, we need about 80 % of the Ay. The ratio A'^"*V^v amounts to about 
0.88 if we employ the value e = 0.01. 



4.2 Ns=32 and 36 Results 

Here we concentrate ourselves to the zero-momentum states which are even under ro- 
tations and reflections. Based on these conditions we construct a basis which possesses 
translation, rotation and reflection symmetries of the lattice [I3]. The total numbers of 
the = basis states are then ~ 2.4 x 10^ and ~ 3.2 x 10^ for the A^s = 32 and 36 
lattices, respectively. Since we can employ any initial trial state as far as it has some 
overlap with the exact eigenstate, we start from the Neel state for simplicity. The value 
of / in the operator Q = II ~ Hj^j,^ is set to be 4.8 (5.4) for the 32-site (36-site) lattice so 
that the largest eigenvalue for Q gives the lowest eigenvalue of Hj^^j^- 

Let us flrst report results on the 32-site lattice. Changing values of e we calculate 
^ j^gg^ fiom 30 samplcs for each value of e. Figs. [7] and [8] present the 
results, where statistical errors are so small that we omit them. In Fig. [7] the results are 
shown as a function of m. We see that the data become almost constant when m >~ 50. 
Fig. E plots values of ((C(io°))) as a function of e. A linear decrease to meet a kink at 
e ~ 0.01 is observed in the Figure, and below this value the results are almost constant. 
We therefore come to a conclusion that the SSS equilibrium manifests itself when e < 0.01 
and that measurements of the normalization factors C*^™-* in this parameter range enable 
us to estimate the lowest energy eigenvalue of the system [22]. We observe that the 
estimated values from ((C(ioo))) with e = 0.01 and e = 0.0025 are E = -16.0052 ± 0.0092 
and E = —16.0047 ± 0.0028, respectively. These values are in good agreement with 
E = —16.0031, the exact ground-state energy obtained from ref. p^. The upper bound 
ofNl""^ is ~ 1.7x10^ for e = 0.01, ~ 1.8xl06for e = 0.005 and ~ l.gxlO^for e = 0.0025, 
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that is about 73%, 77% and 81% of A^y 

Results on the 36-site lattice are qualitatively similar to those on the 32-site lattice. In 
Fig.Elwe show ((^(^o))) for 0.0015 < e < 0.006 as well as {{C'-^^^^)) for 0.0015 < e < 0.0025. 
We observe a linear dependency of {{C^^'^^)) on e in the range 0.003 < e < 0.006, while 
data with less values of e are almost constant. This kind of stability is observed with 
both ((C(40))) and ((C(ioo)^^ fo^. g < 0.0025, although the value of the constant is slightly 
different as we see in the Figure. From 20 samples of ((C(ioo))) ^j^h e = 0.0025, 0.002 
and 0.0015 we obtain E = -18.1317 ± 0.0089, -18.1357 ± 0.0048 and -18.1364 ± 0.0028, 
respectively. Within the statistical errors they all agree with the known exact ground- 
state energy E = —18.1372 [13]. The upper bound of N^^^ (ratio to the A^) is about 
2.4 X 10^ (0.75) when e = 0.0015. 

4.3 Ns = 40 Results 

In this subsection we show results on a 40-site lattice. Since this lattice lacks the reflection 
symmetry, we employ a basis which possesses the up-down symmetry of spins in addition 
to the translation and rotation symmetries of the lattice. The size of the full vector space 
with = is then — 4.3 x 10^. We define Q = 61 — Hj^j^ in this subsection. 

For e > 0.001 we calculate C*-'"-' starting from the Neel state. In the range e > 0.002 
two samples are calculated for each value of e. We observe that deviations of C^"^^ between 
different random number sequences are very small; the statistical errors are 0.5% at most. 
Since our data from one sample would have enough precision for this system, only one 
sample is calculated for each value of e which is less than 0.002. Fig. [10] shows some of 
the results for e > 0.001, where we present data from each sample as a function of m. We 
see that C^"^^ (m > 30) for each value of e becomes almost constant when e > 0.0012. For 
e = 0.001, on the other hand, C*^'"^ is still decreasing at m = 40. We therefore continue 
to calculate C^*") with e = 0.001 up to m = 65. The values we obtain are C^^^^ = 26.620 
and C(65) = 26.417. 

In further calculations of C^"^^ with e = 0.0008, we employ a better trial state which 
is constructed in the same manner as that in ref. [9]. Values of C*-™^ thus obtained are 
presented in Fig. [H] The upper bound of iV^™^ is ~ 3.3 x 10* (76% of A^'v). 

In Fig. [Hwe plot values of C(^°) (for 0.001 < e < 0.0025), the value of C^^^^ (for 
e = 0.001) and the value of C^^o) (for e = 0.0008) function of e. We see that the 
datum with e = 0.0008 is located above the line in the Figure which is determined by 
data between e = 0.001 and e = 0.0015. Therefore, for the same reason as in the 32-site 
and 36-site cases, we conclude that the 40-site system realizes the SSS equilibrium when 
e = 0.0008. 

In order to evaluate the lowest energy of the system without any statistical average, 
we attempt to find a zone which should include the true energy eigenvalue. The upper 
(lower) value of this zone is determined by the maximum (minimum) value of C^'^^ with 
m > rrie, where me denotes a value of m above which the system is in the SSS equilibrium. 
This approach seems to work nicely in the 32-site and the 36-site systems, resultant zones 
being much wider than the stripes obtained by the statistical treatments ^23j. Our best 
estimate for the zone in the 40-site system is —19.92 < E < —19.89, which we obtain 
from values of C^™^ with e = 0.0008 and rrie = 4:0 < m < 50. 
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5 Summary and Discussions 



In this paper we study a frustrated quantum spin model using the SSS (Stochastic State 
Selection) method. Our purpose is to see whether the SSS equilibrium [8], a kind of 
equilibrium which we have found in a study of quantum spin models with positive-definite 
Hamiltonians, is also realized in frustrated systems. 

In the SSS method we start from an initial trial state and recursively calculate the 
m-th normalized intermediate state (m = 1,2,3,---). Each generating procedure is as 
follows. First we operate the m-th random choice operator to the (m — l)-th normalized 
intermediate state. By this operation the effective size of the vector space is drastically 
reduced. The rate of the reduction is controlled by a parameter e. Then successively, we 
operate Q = II — H to the reduced state. This operation again increases the number of 
the basis states which are relevant to the resultant state. Finally we normalize the state 
to obtain the m-th normalization factor and the m-th normalized intermediate state. The 
system is in the SSS equilibrium if, after repeating the procedure many times, the m-th 
normalized intermediate state comes to contain a finite portion of the ground state and 
this portion is irrelevant to m. It should be kept in mind that in unfrustrated models the 
SSS equilibrium is observed for any value of the parameter e. 

What we assert for frustrated quantum spin models is that the SSS equilibrium is also 
realized in these models with small values of the parameter e. After analytical arguments 
in section [3] we present results from numerical calculations for the J1-J2 Heisenberg model 
on 20-, 32-, 36- and 40-site lattices with J2/J1 = 0.5. We observe that these systems are 
in the SSS equilibrium in some range of e. We also confirm that in the equilibrium we 
can effectively estimate the lowest energy eigenvalue of the system. Our results on 20-, 
32- and 36-site lattices are in good agreement with known exact ground-state energies. 
On a 40-site lattice we obtain the result —19.92 < E < —19.89, the basis states required 
for this calculation being maximally 76% of the full vector space. The energy per site 
is then —0.4980 < E/Ng < —0.4973, which is greater than reported values for smaller 
lattices [I3], E/N, = -0.503810, -0.500096 and -0.500615 for A^^ = 36, 32 and 20, 
respectively. Since the value for Ng = 36 seems to be irregular, we omit it from the 
extrapolation. The extrapolation to Ag 00 using our Ag = 20, 32 and 40 data yields 
-0.4993 < (E/N^)^ < -0.4946. 

Several remarks are now in order. First, let us emphasize the wide availability of 
the SSS method in numerical studies. As we have shown in previous works, the SSS 
method is applicable to various models with various symmetries and various system sizes, 
in frustrated cases as well as in unfrustrated ones. 

Using this method it is possible to calculate many physical quantities which are given by 
expectation values of some operators. Essential points of the calculation are as follows. 
Suppose we want to calculate an expectation value {ip \ O \ ip), where O is a given 
operator and | ip) denotes either the exact eigenstate | Ve) or an approximate state 
I iPa). In the SSS method we can generate a state | '?/')+ I x)g instead of | ip), where 
I x)g is the stochastically determined random state. Since | ^/')+ | x)9 is described by 
less number of non-zero coefficients compared to | ip) itself, the method enables us to 
calculate ( {ip \ +g{x | ) O { \ ip)+ \ x')g' )• Then all we should do in addition is to take 
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a statistical average of this quantity, because a mathematical relation 

(( ( I +g{x I ) O ( I ^)+ I ) )) = I O I ^) 

follows from basic properties of the random state. In ref . [7] , for example, we calculate the 
static structure factor of the spin-1/2 anti-ferromagnetic Heisenberg model on a 36-site 
triangular lattice in addition to = and Sz = 1 energy eigenvalues. 

The SSS method provides a fundamental improvement in Monte Carlo techniques. 
Therefore it is easily combined with lots of established techniques such as the power 
method [SIEIC] and the Lanczos method |niin|. A study of the stochastic diagonalization 
by H. De Raedt and M. Frick [21] suggests us that a combination of the SSS method with 
the Jacobi method would also be interesting. 

How large is the cluster size we can manage in calculations with the SSS method, 
then. It of course depends on aims and models. In study of the energy eigenvalue, for 
instance, our purpose in refs. O El [TJ [9] is to obtain upper bounds of eigenvalues from 
approximate eigenstates. While in ref. [8] and in this study we want to evaluate the 
eigenvalue itself through non-zero overlaps with the exact eigenstate which are realized in 
the SSS equilibrium. The system sizes of the latter case are generally smaller than those 
of the former one because smaller values of the parameter e become necessary. 

As for models, in addition to the J1-J2 model we have studied the spin-1/2 Heisenberg 
models on a square lattice up to 64 sites [S], on a triangular lattice up to 48 sites [H], 
the Shastry-Sutherland model up to 64 sites [6]. We found that the method is especially 
effective for the Shastry-Sutherland model. This can be understood because of the very 
compactness of its low-lying eigenstates, where basis states with relatively large coefficients 
dominate if a restructured dimer-like basis |25] is employed. 

We conclude this article with a comment whether one can apply our method to other 
frustrated models. From theoretical points of view, the SSS equilibrium would be observed 
in most models because our discussion in section [3] is based on very moderate assumptions. 
How small the parameter e should be is, however, a model-dependent problem which 
remains to be clarified in future works. If a system realizes the SSS equilibrium with 
moderate values of e it is possible to study many physical properties of the system on 
large lattices. 

Note added 

After submitting this paper for publication we became aware of the work by means of 
the exact diagonalization which suggests that our error estimate might need more 
examinations in future studies. 
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Figure 1: Statistical averages for w^"^^ on the 20-site lattice with e = 0.05 (filled circles), 0.025 (filled 
squares), 0.02 (open squares), 0.015 (open circles) and 0.005 (open diamonds). 
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Figure 2: Statistical averages for C^"*) on the 20-sitc lattice with e = 0.05 (filled circles), 0.03 (filled 
diamonds), 0.025 (filled squares), 0.02 (open squares) and 0.005 (open diamonds). The dashed line 
indicates the eigenvalue Q{= 3 — E) = 13.0123. 
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Figure 3: Statistical averages for (C^™) | | C^™^) (fiUed marks) and (x^^^ | | x^"^) (open marks) 
on the 20-site lattice. Circles, squares and diamonds show the data with e = 0.05, 0.02 and 0.005, 
respectively. 



1.2 



1 X X X X X X X X X X X X X X X X X X X X 



0.8 



0.6 



0.4 



0.2 



-0.2 



qOOOOOOOOOOOOOOOOOO 



o o o o 



50 



100 
m 



150 



200 



Figure 4: Ratios [w^"''>f / [0^""+^^]^ (squares), O2C (l - [w^™^]^) / [C^^+^^j^ (diamonds), 

'32x[fl('")]V[^^''"+'^]' (circles) and {q^ [^i;^]' + Q2C (l - M^^]') + Q2x [s^™)]"} / [C^'^+i)]' 
(crosses) on the 20-site lattice with e = 0.025. 
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Figure 5: Ratios [y,(m)j2 ^ [^(m+i)j2 (squares), Q2C (l - h^™^]^) / [C^™+^^]^ (diamonds), 

Q2x [5'"^]' / (circles) and [q^ [t«M]' + Q2C (l - [w^^^^) + Q^^ [5^"^]'} / [Cf^+i)]' 

(crosses) on the 20-site lattice with e = 0.02. 
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Figure 6: Statistical averages for [5^™^]^ on the 20-site lattice as a function of m. Squares, circles, 
triangles and diamonds show the data with e = 0.02, 0.015, 0.01 and 0.005, respectively. 
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Figure 7: Statistical averages for C^"^^ on the 32-site lattice as a function of m, where each average is 
calculated from 30 samples. Filled diamonds, filled squares, open circles, open squares and open diamonds 
show the data with e = 0.015, 0.013, 0.010, 0.005 and 0.0025, respectively. The dashed line indicates the 
value A.8- E = 20.8031 [13]. 
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Figure 8: Statistical averages for C'-^^"-' on the 32-site lattice as a function of e. Each average is calculated 
from 30 samples. The dashed line indicates the value 4.8 — E — 20.8031 [13 and the dot-dashed line 
shows the result from the linear fit for data with 0.01 < e < 0.015. 
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Figure 9: Statistical averages for (7^'*°-' (open squares, open diamonds) and (7^^°"' (filled diamonds) on the 
36-site lattice as a function of e. Each average for e < 0.0025 (the open diamond or the filled diamond) 
is calculated from 20 samples, while each averages for e > 0.003 (the open square) is calculated from 2 
samples . The dashed line indicates the value of 5.4 — E — 23.5372 [13^. The dot-dashed line shows the 
result from the least square fit for data with 0.003 < e < 0.006. 
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Figure 10: Values of C^'") on the 40-site lattice for e > 0.001 as a function of m. The initial trial state is 
the Neel state. Open circles, filled circles, open squares, filled squares, open diamonds, filled diamonds 
and asterisks denote data with e = 0.005, 0.0035, 0.0025, 0.002, 0.0015, 0.0012 and 0.001, respectively. 
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Figure 11: Values of C^"*' on the 40-site lattice with e = 0.0008 as a function of m. The initial trial state 
is an approximate ground state which yields (^-^"^ | Q \ tp^°^) = \6i -H \ tp^°^) = 18.814. 
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Fiffurc 12: Values of C^^o) (open triangles), C^'™' (the filled square) and C^^s) (the open circle) on the 
40-site lattice as a function of e. Each datum is obtained from one sample. The dot-dashed line shows 
the result from the least square fit for data with 0.001 < e < 0.0015. 
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